Simplex solids in SU(iV) Heisenberg models on the kagome and checkerboard lattices 
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We present a numerical study of the SU(iV) Heisenberg model with the fundamental representation at each 
site for the kagome lattice (for N = 3) and the checkerboard lattice (for TV = 4), which are the line graphs of 
the honeycomb and square lattices and thus belong to the class of bisimplex lattices. Using infinite projected 
entangled-pair states (iPEPS) and exact diagonalizations, we show that in both cases the ground state is a simplex 
solid state with a two-fold degeneracy, in which the N spins belonging to a simplex (i.e. a complete graph) form 
a singlet. Theses states can be seen as generalizations of valence bond solid states known to be stabilized in 
certain SU(2) spin models. 
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I. INTRODUCTION 

The SU(2) Heisenberg model on the square lattice is one 
of the most studied systems in condensed matter physics, and 
its properties are by now well understood. Generalizations of 
this model to SU(iV) with different values of N, lattice ge- 
ometries, and representations of SU(iV) exhibit an extremely 
rich variety of different ground states. While such models 
have been the subject of many theoretical studies in the past 
decades, 1-9 they have recently attracted increasing interest 
thanks to the proposals to realize SU(iV) symmetric Hubbard 
models in experiments on ultracold fermionic atoms in optical 
lattices. 10 " 13 

As in the SU(2) case, the SU(N) Heisenberg models can 
be derived from a generalized Hubbard model of N flavors 
of fermions in a Mott insulating phase. Different numbers 
of particles per site correspond to different representations of 
SU(iV). In this work we focus on the Heisenberg model with 
one particle per site, corresponding to the fundamental rep- 
resentation (a Young diagram with a single box) of SU(iV) 
at each site. The Hamiltonian can be written in terms of a 
permutation operator Pij 14 which exchanges the particles on 
neighboring sites, 



(i) 



(id) 



The expectation value of the operator Py is minimal (and 
equal to —1) when the wave function is antisymmetric on the 
ij bond. For SU(iV), a fully antisymmetric wave function can 
be constructed using exactly N sites (and not only two sites as 
in the familiar SU(2) case). This is an SU(A r ) singlet, and the 
energy on any bond belonging to the singlet reaches its lowest 
possible value — 1 . This should be contrasted to the other class 
of SU{N) Heisenberg models, 2 7 where conjugate representa- 
tions on two different sublattices were used, and thus singlets 
can be formed between two sites for any 7Y. 
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FIG. 1 : (Color online) The simplex solid states obtained with iPEPS 
(D = 14) for two different SU(7V) Heisenberg models. The width 
of a bond is proportional to the magnitude of the bond energy, 
while blue (red dotted) bonds correspond to negative (positive) en- 
ergies, (a) One of the two degenerate trimerized ground states of the 
SU(3) Heisenberg model on the kagome lattice, (b) One of the two 
quadrumerized ground states of the SU(4) Heisenberg model on the 
checkerboard lattice. 



In this work we study the model on a particular class of 
lattices, so-called bisimplex lattices, 15 consisting of corner 
sharing simplices residing on an underlying bipartite lattice 
(equivalent to the line graph of the underlying bipartite lat- 
tice). In particular, we investigate the possibility of having 
an 7V-merized ground state, i.e. where the lattice is cov- 
ered by singlets which extend over TV-site simplices. 48 Such 
simplex solid states 49 can be seen as a generalization of the 
valence bond solids (or valence bond crystals) known from 
certain SU(2) models, as for example the dimerized state in 
the Majumdar-Gosh model. Examples we consider here are 
the kagome lattice with 3-site simplices (triangles), and the 
checkerboard lattice with 4-site simplices (crossed-squares in 
Fig. 1(b), which are topologically equivalent to tetrahedra). 
Because of the underlying bipartite lattice (honeycomb for 
kagome, square for checkerboard) there are two possible ways 
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(a) Setup 1 (b) Setup 2 




FIG. 2: (Color online) The two different simulation setups used to 
simulate the kagome lattice with the iPEPS method developed for 
the square lattice. Black circles and triangles correspond to tensors, 
black lines to connection between tensors (the physical index of a 
tensor is not shown). Interactions between physical sites (filled cir- 
cles) are given by thick shaded lines, (a) Auxiliary tensors (white cir- 
cles) are introduced to create a square lattice iPEPS. The interactions 
along the horizontal and vertical direction correspond to nearest- 
neighbor couplings on the square lattice, whereas the interactions 
along the diagonal are treated as next-nearest neighbor interactions, 
(b) Three physical sites A, B, C on the kagome lattice, each having 
a local dimension d — 3, are mapped into a block site with a local 
dimension d — 27. 



to cover the lattice with singlets leading to a two-fold ground 
state degeneracy. 

The occurrence of a symmetry broken state with a two- 
fold degeneracy on these lattices has already been observed 
in Refs. 16,17 for a t-J model away from half filling, which 
showed that this type of symmetry breaking is a likely can- 
didate on these lattices. In Ref. 18 Arovas derived parent 
Hamiltonians for the exact SU(A r ) simplex solid states on the 
kagome (for N = 3) and the checkerboard lattice (for N = 4). 
It is conceivable that these parent Hamiltonians are adiabati- 
cally connected to the SU( AO Heisenberg model, in the same 
way as the Affleck-Kennedy-Lieb-Tasaki (AKLT) state 19 is 
adiabatically connected with the ground state of the 5=1 
spin chain. Simplex solid states on bisimplex lattices have 
also been predicted in Ref. 9 based on studies of a particular 
large N limit, with representations labelled by Young tableaus 
with m rows and n c columns, where N/m = k and n c are 
held fixed. 

We use infinite projected entangled-pair states (iPEPS) and 
exact diagonalization (ED) to study the ground state of the 
N = 3 case on the kagome and the N = 4 case on the 
checkerboard lattice. Both methods consistently predict a 
two-fold degenerate simplex solid ground state for both mod- 
els, summarized in Fig. 1. 

The paper is organized as follows: In Sec. II we provide de- 
tails on the iPEPS simulations, in particular, how the models 
are simulated using a square-lattice iPEPS. In Sees. Ill and IV 
we present the iPEPS and ED results obtained for the kagome 
and the checkerboard model, respectively. Finally, Sec. V 
summarizes our findings. 



H. INFINITE PROJECTED ENTANGLED-PAIR STATES 
(IPEPS) 

A projected-entangled pair state (PEPS) is an efficient vari- 
ational ansatz for two-dimensional ground state wave func- 
tions. 20-24 It can be seen as a natural extension of a matrix 
product state (MPS), the underlying ansatz of the famous den- 
sity matrix renormalization group (DMRG) method. 25 The 
main idea is to represent a wave function by a trace of a prod- 
uct of tensors, with one tensor per lattice site. On the square 
lattice each tensor Tf- kl has five indices: one index p which 
carries the local Hilbert space of a lattice site with dimension 
d, and four indices i, j, k, I - the auxiliary bonds with bond di- 
mension D - which connect to the four nearest-neighbor ten- 
sors. Thus, each tensor contains dD 4 variational parameters, 
and by varying D the accuracy of the ansatz can be system- 
atically controlled. A bond dimension D = 1 simply cor- 
responds to a product state (a site-factorized wave function), 
and upon increasing D quantum fluctuations (entanglement) 
can be taken into account in a systematic way. 

Details on the iPEPS method for the square lattice can be 
found in Refs. 26-28, in particular how to optimize the ten- 
sors (i.e. finding the best variational parameters) and how 
to compute expectation values by contracting the tensor net- 
work (i.e. computing the trace of the product of all ten- 
sors). We performed similar iPEPS simulations already for 
the SU(4) Heisenberg model 29 and the triangular- and square- 
lattice SU(3) Heisenberg model. 30 

For the experts, we note that the optimization is done 
through an imaginary time evolution using the simple 
update, 27 ' 3 1-33 and we have verified some simulation re- 
sults also with the full update. 27 The corner-transfer matrix 
method 27 34 ' 35 is used to contract the tensor network, where 
the error of the approximate contraction can be controlled by 
the boundary dimension \- T ne simulation results in this work 
are extrapolated in where the extrapolation error is small 
compared to the symbol sizes. To improve the efficiency of 
the simulations we used tensors with 7L q symmetry. 36,37 

To simulate the checkerboard model we use the usual 
square lattice iPEPS ansatz where we treat the diagonal cou- 
plings as next-nearest neighbor interactions as described in 
Ref. 38. 

For the kagome lattice we implemented two different simu- 
lation setups based on a square lattice iPEPS, which have the 
advantage that existing algorithms for the optimization and 
contraction can be reused. For the first variant we use one 
tensor per lattice site, plus additional auxiliary tensors which 
are inserted to form a square lattice iPEPS, as sketched in 
Fig. 2(a). [The bond dimension of the auxiliary tensors can 
be chosen as D = 1 since all correlations are carried by the 
tensors on the physical sites.] The couplings along the hor- 
izontal and vertical direction correspond to nearest-neighbor 
couplings between two tensors, whereas the remaining bonds 
necessary to form the kagome lattice are represented by the 
next-nearest-neighbor bonds in the square lattice, which can 
be treated as explained in Ref. 38. In the second setup we map 
the kagome lattice onto a square lattice by blocking three sites 
as illustrated in Fig. 2(b). The original Hamiltonian is mapped 



3 



-0.2 
-0.4 
lu™ -0.6 
-0.8 
-1 

0.8 
0.6 
E 0.4 
0.2 





(a) 


□ iPEPS setup 1 






— e— iPEPS setup 2 






* ED 











0.5 




0.1 0.2 0.3 0.4 
1/D OI-1/N 



0.1 0.2 0.3 0.4 

1/D 



(C) 





0.1 0.2 0.3 0.4 

1/D 



FIG. 3: (Color online) (a-c) iPEPS simulation results for the SU(3) 
Heisenberg model on the kagome lattice as a function of inverse bond 
dimension 1/D for the two simulation setups (cf. Fig. 2). Extrap- 
olations (dashed lines) are a guide to the eye. (a) Energy per site 
compared with the results from ED (plotted as a function of 1/N S ). 
The estimated value in the infinite D limit is —0.829(1). (b) The 
difference in bond energies AE remains finite in the infinite D limit, 
which shows that the state breaks translational invariance as illus- 
trated in Fig. 1(a). (c) Local moment which is strongly suppressed 
with increasing bond dimension, (d) ED results for the energy gaps 
to the first singlet and the first SU(3) multiplet excitations plotted as 
a function of the inverse number of sites 1/N 3 . 



onto a square lattice Hamiltonian with nearest-neighbor inter- 
actions between the block sites (see supplementary material). 
We point out here that we do not block three sites belonging to 
a triangle in the kagome lattice, since this would automatically 
bias the solution towards a trimerized state. 

Since we work directly in the thermodynamic limit, the 
ground state wave function may exhibit spontaneously bro- 
ken symmetries. We measure the energy on each bond, Ej,, 
in the unit cell. If the energies are not equal on all symmetry 
related bonds, i.e. if the difference 



AE = max(.E&) — mm(Eb) 



(2) 



is finite, then the state breaks some lattice symmetries. 

Furthermore, we compute the local ordered moment m on 
each site, 



m 



N 



N 



1 ^ 



(3) 



where S@ = |a)(/3| are the generators of SU(A) and a, f3 
run over all flavor indices. A finite m implies that the SU(A) 
symmetry is broken. 

We tried different unit cell sizes up to size 6 x 6 in iPEPS 
to test for possible stable color-ordered solutions, however all 
simulations consistently lead to an A-merized ground state. 
The results presented in the following are obtained with a 2 x 2 
unit cell of tensors. 



III. THE SU(3) MODEL ON THE KAGOME LATTICE 

We consider the model ( 1 ) on the kagome lattice for N = 3 
with the fundamental representation at each lattice site. Sin- 
glet formation most naturally occurs between three sites on a 
triangle, resulting into two possible coverings on the kagome 
lattice in which either all triangles pointing upwards, or all 
triangles pointing downwards form singlets. The latter case is 
illustrated in Fig. 1(a). 

Figure 3(a) shows the iPEPS energy per site, E s obtained 
with the two different simulation setups (cf. Fig. 2) as a func- 
tion of bond dimension D. It is in general not known how 
quantities depend on D. Here we find that for D > 10 the en- 
ergy does not change much anymore upon further increasing 
D. As an estimate in the infinite D limit we take the middle 
between the value at the maximal D and the value obtained 
from a linear extrapolation in 1/D for the four largest D. 

From Fig. 3(b), it is clear that the difference between the 
highest and lowest bond energy AE is finite for all values of 
D, which indicates that the ground state breaks translational 
symmetry. The strength of the individual bond energies in 
the lattice can be seen in Fig. 1(a), where the thickness of a 
bond is proportional to the magnitude of the bond energy. The 
resulting pattern is clearly compatible with a trimerized state. 
The strength of the trimerization decreases with increasing D 
(except for the lowest D = 2), but it seems to remain finite 
around AE = 0.76(1) in the infinite D limit. 

Since the trimerized state does not break the SU(3) symme- 
try, the local ordered moment m defined in Eq. (3) should van- 
ish. Figure 3(c) shows that m is strongly suppressed with in- 
creasing D. In both simulation setups m vanishes completely 
for D > 10. 

We note that we have also tested a 6 x 6 unit cell which 
can accommodate the \/3 x \/3 color ordered state predicted 
by linear flavor-wave theory (see Sec. V). We found that this 
state only appears as a metastable state in iPEPS, i.e., the N- 
merized ground state has a lower variational energy. 

To further corroborate the iPEPS results we have performed 
ED simulations of the SU(3) Heisenberg model on finite size 
kagome samples of up to N s = 24 sites (cf. Ref. 39 re- 
garding the geometry of the clusters). The energies per site 
have been included in Fig. 3(a) and the agreement between 
the two methods at large D or N s is very good. In Fig. 3(d) 
we display two relevant energy gaps, the first one to the low- 
est singlet excitation at the T point, and the second one to 
the first magnetic (i.e. non-singlet) excitation. While the sin- 
glet gap seems to collapse for larger N s , in agreement with 
the lattice symmetry breaking scenario requiring a two-fold 
ground state degeneracy, the magnetic gap appears to stay fi- 
nite, in agreement with the picture of a nonmagnetic singlet 
ground state. In the left panel of Fig. 4 we finally show the 
connected bond energy correlations (cf. caption), which con- 
vincingly demonstrate the trimerization of the kagome lattice 
SU(3) Heisenberg model. 

In summary, both iPEPS and ED provide strong indications 
that the ground state of the SU(3) Heisenberg model on the 
kagome lattice is trimerized. 
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FIG. 4: (Color online) ED results for the connected bond energy cor- 
relations (PijPki) — (Pij) {Phi} - The reference bond is black, pos- 
itive (negative) correlations are shown with blue (red dashed) lines. 
The width of the lines is proportional to the correlation function. The 
periodic (Wigner-Seitz) cell is shaded in grey. Left panel: Results for 
the SU(3) Heisenberg model for a N s = 21 kagome sample. Right 
panel: Results for the SU(4) Heisenberg model on the checkerboard 
lattice for a N s = 20 sample. 



IV. THE SU(4) MODEL ON THE CHECKERBOARD 
LATTICE 

The SU(4) model on the checkerboard lattice is described 
by the Hamiltonian ( 1 ) with the fundamental representation of 
SU(4) at each lattice site. Thus the local dimension of a lattice 
site is four. 

Figure 5(a) shows the iPEPS energies as a function of in- 
verse bond dimension. We find that even for the largest values 
of D used the energy still decreases considerably upon further 
increasing D. In the infinite D limit we expect the energy to 
lie in the range -1.45 < E a < -1.30. The ED result for 
finite systems E**= 16 = -1.392, and £f== 20 = -1.347 lie 
within this range. 

The difference in bond energies AE, shown in Fig. 5(b), 
becomes large for D > 5. Not all the weak bonds have 
exactly the same energy, see Fig. 5(d). So we took aver- 
aged energies of the weak and strong bonds to compute AE. 
The distribution of the individual bond energies on the lat- 
tice can be seen in Fig. 1(b), clearly showing that quadrumers 
are formed. With increasing bond dimension AE slightly de- 
creases but the data still suggest a rather large value in the 
range AE « 0.77 - 0.80 in the infinite D limit. 

Figure 5(c) shows that the ordered moment vanishes in the 
large D limit, compatible with a quadrumerized state which 
does not break SU(4) symmetry. 

Finally, ED results for the connected bond energy correla- 
tions are shown in the right panel of Fig. 4 for a N s = 20 sam- 
ple. They further highlight the strong quadrumerization insta- 
bility of the checkerboard lattice SU(4) Heisenberg model. 



V. DISCUSSION 

The iPEPS and exact diagonalization results reported in this 
paper give clear evidence in favor of the stabilization of sim- 
plex solids for SU(3) on kagome and for SU(4) on the checker- 
board lattice. The underlying mechanism is the formation of 
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FIG. 5: (Color online) iPEPS simulation results for the SU(4) 
Heisenberg model on the checkerboard lattice as a function of in- 
verse bond dimension, (a) Energy per site, E 3 , compared with the 
results from ED (plotted as a function of 1/N 3 ). (b) Different bond 
energies in the quadrumerized state. The strong horizontal and verti- 
cal bonds (H/V) have the same energy as the strong diagonal bonds 
for large D. (c) The difference in bond energies AE is finite which 
shows that the state breaks lattice symmetries [cf. Fig. 1(b)]. (d) The 
local moment decreases rapidly with D and vanishes for large D. 



SU(N) singlets on the simplices of the lattice, which is pos- 
sible whenever N is equal to the number of sites of the sim- 
plices. Based on these results we also expect a two-fold de- 
generate simplex solid ground state for the SU(4) model on 
the pyrochlore lattice (as suggested by previous work 15 ' 16 ' 18 ) 
as well as the SU(3) model on the Garnet lattice or variants 
thereof. 15 - 40 

In principle, it is sufficient that the number of sites be a 
multiple of N for this scenario to be realized. Interestingly, 
the SU(2) Heisenberg model on the checkerboard lattice ful- 
fills this condition, and its ground state is also spontaneaously 
quadrumerized. 41 ' 42 However, the strong bonds form plaque- 
ttes on the empty squares in that case, and the mechanism is 
different: the diagonal bonds are instrumental in frustrating 
the inter-plaquette coupling. 

Another interesting conclusion can be drawn regarding 
flavor-wave theory. In models investigated so far, flavor-wave 
theory (FWT) has always revealed some aspect of the ground 
state: it led to the correct ground state for the SU(3) model on 
the square 43 and the triangular 44 lattices, and it gave interest- 
ing insight for the SU(4) model on the square lattice. 29 In the 
case of SU(3) on kagome reported in this paper, it fails com- 
pletely: the ground state is predicted to be an ordered state 
with a \/3 x \/3 unit cell and no trimerization (all bond en- 
ergies remain equal). In fact, it is simply impossible to de- 
scribe a trimerized state with FWT since this would required 
to build an unfrustrated, two-color configuration on a triangle. 
It is nevertheless amusing to notice that the ground state de- 
generacy is already lifted at the harmonic level, by contrast to 
the SU(2) case, and that the \/3 x \/3 state (energy per site of 
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— 2 + 2/^/3 ~ —0.845) is favored over the q = state (energy 
per site of — 2 + 4/ir ~ —0.727), as for the SU(2) case beyond 
linear spin-wave theory. 45 ^ 7 We note that these energies are 
not variational and thus cannot be compared with the iPEPS 
energies. 
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I. DETAILS ON THE SETUP 2 TO SIMULATE THE 
KAGOME MODEL 

Here we explain how to map a nearest-neighbor Hamilto- 
nian on the kagome lattice h onto a nearest-neighbor Hamil- 
tonian on the square lattice H, where three physical sites A, 
B, C on the kagome lattice are blocked into one block site of 
the square lattice, as shown in Fig. 2(b) in the main text. The 
Hilbert space of a block site Hbs is given by 



Ubs =Ua®Ub®'Hc, 



(1) 



where Ha, T~Lb, He are the local Hilbert spaces of the three 
physical sites A, B, C, respectively. Thus, for the SU(3) 
Heisenberg model with local physical dimension 3 the local 
dimension of a block site is 3 3 = 27. 

The original Hamiltonian defined on the kagome lattice is 
given by 



(2) 



where the sum (u, v) goes over nearest-neighbor sites on the 
kagome lattice. The Hamiltonian H defined on the lattice 
formed by the block sites reads 



(3) 



where the sum i goes over all block sites in the square lattice, 
and (i,j) x and (i,j) y are sums over nearest-neighbor pairs in 
x-direction and y-direction, respectively. Each site term Hf 
contains two terms of the original Hamiltonian (2): 



Hf = h t 



(4) 



where hi 



(hi 



,) is the Hamiltonian term connecting the 



physical sites A and B ( A and C) on block site i. These 



two terms correspond to the horizontal and vertical blue lines 
within a block site in Fig. 2(b). In x-direction the block 
Hamiltonian is given by the sum of two contributions, 



TJX 



(5) 



where is is the i?-site on block site i, and ja is the A-site 
of the neighboring block site j. The first (second) term corre- 
sponds to the horizontal (diagonal) blue line in Fig. 2(b) which 
connect two neighboring block sites in a>direction. Similarly, 
the block Hamiltonian term in y-direction reads 



H y - h 

i,j ~ l c,3A 



^CJB 



(6) 



where the first (second) term corresponds to the vertical (di- 
agonal) blue line in Fig. 2(b) which connect two neighboring 
block sites in y-direction. 

Thus, we find that the Hamiltonian (3) defined on the block 
site lattice is a nearest-neighbor Hamiltonian on the square lat- 
tice, which can be simulated with the standard iPEPS method 
for the square lattice. 

We point out that the correlations within a block site are 
taken into account in an exact way, whereas the correlations 
between physical sites on neighboring block sites depend on 
the bond dimension D. Thus, by performing the blocking we 
favor correlations between the three physical sites in a block 
site, which can bias the result for small bond dimensions. 
However, as we note in the paper, we do not block the three 
sites which eventually form a singlet (if we would do this we 
would bias the state to a singlet state from the start, i.e. a triv- 
ial D = 1 PEPS could already represent a trimerized state). 
This bias becomes irrelevant for large bond dimensions, be- 
cause for large D the physical sites on neighboring block sites 
can also become strongly correlated. Indeed, as we show in 
the paper, at large bond dimensions the results become equiv- 
alent with the ones obtained with the simulation setup I, where 
we did not use this blocking. 



